Topological phase separation in 2D hard-core Bose-Hubbard system away from 

half-filling 



O 

o 

(N 

< 
(N 



A.S. Moskvin, 1 I.G. Bostrem, 1 and A.S. Ovchinnikov 1 

1 Ural State University, 620083, Ekaterinburg, Russia 
(Dated: February 2, 2008) 

We suppose that the doping of the 2D hard-core boson system away from half-filling may result 
in the formation of multi-center topological defect such as charge order (CO) bubble domain(s) with 
Bose superfluid (BS) and extra bosons both localized in domain wall(s), or a topological CO+BS 
phase separation, rather than an uniform mixed CO+BS supersolid phase. Starting from the classical 
model we predict the properties of the respective quantum system. The long-wavelength behavior of 
the system is believed to remind that of granular superconductors, CDW materials, Wigner crystals, 
and multi-skyrmion system akin in a quantum Hall ferromagnetic state of a 2D electron gas. 
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One of the fundamental hot debated problems in 
bosonic physics concerns the evolution of the charge or- 
dered (CO) ground state of 2D hard-core Bose-Hubbard 
(BH) model (hc-BH) with a doping away from half-filling. 
Numerous model studies steadily confirmed the emer- 
gence of "supersolid" phases with simultaneous diago- 
nal and off-diagonal long range order while Penrose and 
OnsageitA' were the first showing as early as 1956 that 
supersolid phases do not occur. 

The most recent quantum Monte-Carlo (QMC) 
simulations 2,3 '! found two significant features of the 2D 
BH model with a screened Coulomb repulsion: the ab- 
sence of supersolid phase at half-filling, and a grow- 
ing tendency to phase separation (CO+BS) upon dop- 
ing away from half-filling. Moreover, Batrouni and 
Scalettar— studied quantum phase transitions in the 
ground state of the 2D hc-BH Hamiltonian and have 
shown numerically that, contrary to the generally held 
belief, the most commonly discussed "checkerboard" su- 
persolid is thermodynamically unstable and phase sep- 
arates into solid and superfluid phases. The physics of 
the CO+BS phase separation in BH model is associated 
with a rapid increase of the energy of a homogeneous CO 
state with doping away from half-filling due to a large 
" pseudo-spin- flip" energy cost. Hence, it appears to be 
energetically more favorable to "extract" extra bosons 
(holes) from the CO state and arrange them into finite 
clusters with a relatively small number of particles. Such 
a droplet scenario is believed to minimize the long-range 
Coulomb repulsion^ However, immediately there arise 
several questions. Whether a simple mean-field approx- 
imation (MFA) and classical continuum model can pre- 
dict such a behavior? What is the detailed structure 
of the CO+BS phase separated state? What are the 
main factors governing the essential low-energy and long- 
wavelength physics? Is it possible to make use of simple 
toy models? 

In the Letter we present a topological scenario of 
CO+BS phase separation in 2D hc-BH model with inter- 
site repulsion. The extra bosons/holes doped to a 
checkerboard CO phase of 2D boson system are believed 
to be confined in the ring-shaped domain wall of the 
skyrmion-like topological defect which looks like a bub- 



ble domain in an easy- axis antiferromagnet. This allows 
us to propose the mechanism of 2D topological CO+BS 
phase separation when the doping of the bare checker- 
board CO phase results in the formation of a multi-center 
topological defect, which simplified pseudo-spin pattern 
looks like a system of bubble CO domains with Bose su- 
perfluid confined in charged ring-shaped domain walls. 

The Hamiltonian of hard-core Bose gas on a lattice can 
be written in a standard form as follows: 



Hbg = -J2 t v P (BjB J +B}B l )P- 
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where P is the projection operator which removes double 
occupancy of any site, Ni — BjBi, fj, the chemical poten- 
tial determined from the condition of fixed full number of 

N 

bosons Ni — (N) or concentration n = Ni/N £ [0, 1]. 
»=i 

The tij denotes an effective transfer integral, Vij is an 
intersite interaction between the bosons. Here B'(B) 
are the Pauli creation (annihilation) operators which are 
Bose-like commuting for different sites [Bi,B^] = 0, 
for i j, while for the same site Bf — (B\) 2 = 0, 
[Bi, B\] = 1 - 2N l7 Ni = B\Bi, N is a full number of 
sites. 

The model of hard-core Bose-gas with an intersite re- 
pulsion is equivalent to a system of spins s = 1/2 exposed 
to an external magnetic field in the ^-direction 
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where Jf? 



(2) 



It T z 



Vij, h 



= M - E Vij, s = 

^= 2 B,s+ = -^= 2 B\s* = -\+B\Bi, a± = t^(**±m»). 

In frames of of a conventional two-sublattice MFA ap- 
proach one may introduce two vectors, the ferromagnetic 
and antiferromagnetic ones: 

m = + (s 2 »; 1 = f «si) - (s 2 )); m 2 + l 2 = 1 . 

2s Zs 

where m ■ 1 = 0. For the plane where these vec- 
tors lie one may introduce two-parametric angular de- 
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scription: m x = sin a cos /3, m z = — sin a sin /3, = 
cos a sin /3, l z — cos a cos/?. The hard-core boson sys- 
tem in a two-sublattice approximation is described by 
two diagonal order parameters l z , m z and two complex 
off-diagonal order parameters m± = =p-^=(m x ± vm y ) 

and Z± = -F^Qx ± ^y)- The complex superfluid order 
parameter ^(r) = ^(r)! exp — up is determined by the 
in-plane components of ferromagnetic vector: W(r) = 
\{{B\ + B2)) = sm_ = smj_ exp— m± being the 
length of the in-plane component of ferromagnetic vec- 
tor. So, for a local condensate density we get n s — s 2 m 2 L . 
It is of interest to note that in fact all the conventional 
uniform T = states with nonzero \P(r) imply simul- 
taneous long-range order both for modulus |\£(r)| and 
phase tp. The in-plane components of antifcrromagnetic 
vector l± describe a staggered off-diagonal order. It is 
worth noting that by default one usually considers the 
negative sign of the transfer integral ty , that implies the 
ferromagnetic in-plane pseudospin ordering. 

The model exhibits many fascinating quantum phases 
and phase transitions. Early investigations predict at 
T = charge order (CO), Bose superfluid (BS) and 
mixed (CO+BS) supersolid uniform phases with an 
Ising-type melting transition (CO-NO) and Kosterlitz- 
Thouless-type (BS-NO) phase transitions to a non- 
ordered normal fluid (NO) A The detailed mean- field and 
spin-wave analysis of the uniform phases of 2D hc-BH 
model is given by Pich and Freyi 

MFA yields for the conventional uniform supersolid 
phase 7 



. 2a y/V=2t . 2 

sin 5 = m z . ; sin a 

VV + 2t 
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with a constant chemical potential \x — 4s^/ (V 2 — At 2 ). 
It should be noted that the supersolid phase appears to 
be unstable with regard to small perturbations in the 
Hamiltonian. The mean-field energy per site of the uni- 
form supersolid phase is written as follows^ 

Ess = E C o + s/im 2 = E C o + K n B ~\)t 

where Eco = — 2Vs 2 . The cost of doping both for CO 
and CO+BS phase appears to be rather high as com- 
pared with an easy-plane BS phase at half-filling where 
the chemical potential turns into zero£ 

Magnetic analogy allows us to make unambiguous pre- 
dictions as regards the doping of BH system away from 
half-filling. Indeed, the boson/hole doping of checker- 
board CO phase corresponds to the magnetization of an 
antiferromagnet in z-direction. In the uniform easy-axis 
Z z -phase of anisotropic antiferromagnet the local spin-flip 
energy cost is rather big. In other words, the energy cost 
for boson/hole doping into checkerboard CO phase ap- 
pears to be big due to a large contribution of boson-boson 
repulsion. 

However, the magnetization of the anisotropic anti- 
ferromagnet in an easy axis direction may proceed as 



a first order phase transition with a "topological phase 
separation" due to the existence of antiphase domains. 9 
The antiphase domain walls provide the natural nucle- 
ation centers for a spin-flop phase having enhanced mag- 
netic susceptibility as compared with small if any lon- 
gitudinal susceptibility thus providing the advantage of 
the field energy. Namely domain walls would specify 
the inhomogeneous magnetization pattern for such an 
anisotropic easy-axis antiferromagnet in relatively weak 
external magnetic field. As concerns the domain type 
in quasi-two-dimensional antiferromagnet one should em- 
phasize the specific role played by the cylindrical, or bub- 
ble domains which have finite energy and size. These 
topological solitons have the vortex-like in-plane spin 
structure and often named " skyrmions" . The classical, 
or Belavin-Polyakov (BP) skyrmioniS describes the solu- 
tions of a non-linear cr-model with a classical isotropic 2D 
Heisenberg Hamiltonian and represents one of the generic 
toy model spin textures. It is of primary importance to 
note that skyrmion-like bubble domains in easy-axis lay- 
ered antiferromagnets were actually observed in the ex- 
periments performed by Waldner^i that were supported 
later by different authors (see, e.g. Refs ll2lll3lK Al- 
though some questions were not completely clarified and 
remain open until nowji^i^ the classical and quantum^ 
skyrmion-like topological defects are believed to be a gen- 
uine element of essential physics both of ferro- and anti- 
ferromagnetic 2D easy-axis systems. 

The magnetic analogy seems to be a little bit naive, 
however, it catches the essential physics of doping the 
hc-BH system. As regards the checkerboard CO phase of 
such a system, namely a finite energy skyrmion-like bub- 
ble domainii seems to be the most preferable candidate 
for the domain with antiphase domain wall providing the 
natural reservoir for extra bosons. The classical descrip- 
tion of nonlinear excitations in hc-BH model includes the 
skyrmionic solution given V = 2t*L The skyrmion spin 
texture consists of a vortex-like arrangement of the in- 
plane components of ferromagnetic m vector with the l z - 
component reversed in the centre of the skyrmion and 
gradually increasing to match the homogeneous back- 
ground at infinity. The simplest spin distribution within 
classical skyrmion is given as follows 

m x = to j_ cos((y3 + tpo); m y — m± sia(cp + ipq); 



m j_ 



2r\ 
r 2 + A 2 ' 



h = 



\ 2 



(3) 



where ipo is a global phase (U(l) order parameter), A 
skyrmion radius. The skyrmion spin texture describes 
the coexistence and competition of the charge order pa- 
rameter l z and BS order parameter m± that reflects a 
complex spatial interplay of potential and kinetic ener- 
gies. The skyrmion looks like a bubble domain in an 
easy-axis magnet. Skyrmions are characterized by the 
magnitude and sign of its topological charge, by its size 
(radius), and by the global orientation of the spin. The 
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scale invariance of classical BP skyrmionic solution re- 
flects in that its energy does not depend on radius, and 
global phase. An interesting example of topological inho- 
mogeneity is provided by a multi-center BP skyrmioniS 
which energy does not depend on the position of the cen- 
ters. The latter are believed to be addressed as an addi- 
tional degree of freedom, or positional order parameter. 

It should be noted that the domain wall in such a 
bubble domain somehow created in the checkerboard CO 
phase of 2D hc-BH system represents an effective ring- 
shaped reservoir both for Bose-superfluid and extra bo- 
son/hole. Indeed, the magnetic response of a fixed radius 
skyrmion to a weak longitudinal external magnetic field 
can be approximated as follows: 

m z (r) w x( r ) h , 

where x( r ) — /( r ) m J-( r ) is a local magnetic susceptibil- 
ity, and f(r) has distinct maximum at the domain wall 
center r = A with oscillating behavior f(r) cx J\(r\/2) oc 
sin(rv / 2 + ir/A)/y/r for large r. The oscillating charac- 
ter of the magnetization distribution seems to reflect a 
skyrmion instability. Thus we conclude that the extra 
boson/hole density described by the order parameter m z 
is to be localized in the domain wall. This result is sup- 
ported by numerical calculations for the static skyrmion- 
like excitations in a continuous 2D hc-BH model with 
a strong boson-boson repulsion that was considered by 
us earlier in The Fig. 1 sketches the calculated 

radial distribution of the order parameter m z for a sin- 
gle bubble domain in hc-BH system with injected small 
boson concentration away from half-filling. It appears 
to be trapped inside the domain wall. As expected, the 
soliton energy depends quadratically on the number An 
of bosons bound in domain wall (see insert in Fig. 1), 
similarly to that of homogeneous BS phase. 8 In other 
words, dE/dn — and one might say about a zero value 
of the effective boson/hole chemical potential for the CO 
bubble domain configuration would it be a ground state. 

In the continuous model the classical BP skyrmion is 
a topological excitation and cannot dissipate. However, 
the classical static skyrmion is unstable with regard to an 
external field, anisotropic interactions both of easy-plane 
and easy-axis type. Small easy-axis anisotropy or exter- 
nal field are sufficient to shrink skyrmion to a nanoscopic 
size when magnetic length Iq: 

is of the order of several lattice parameters, and the 
continuous approximation fails to correctly describe ex- 
citations. Nonetheless, Abanov and Pokrovskyi 8 . have 
shown that the easy-axis anisotropy together with fourth- 
order exchange term can stabilize skyrmion with radius 

R CX y/lQ. 

A skyrmionic scenario in hc-BH model allows us to 
make several important predictions. Away from half- 
filling one may anticipate the nucleation of a topologi- 
cal defect in the unconventional form of the multi-center 
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FIG. 1: The results of numerical calculations by shooting 
method for a shrinked skyrmion given v = 2.1, t = 1: main 
panel - the radial distribution of the order parameter m z (the 
filled region points to a center of a domain wall); the insert 
shows the dependence of the soliton energy on the number 
An of bosons bound in a domain wall. 



skyrmion-like object with ring-shaped Bose superfluid re- 
gions positioned in an antiphase domain wall separating 
the CO core and CO outside of the single skyrmion. The 
specific spatial separation of BS and CO order parame- 
ters that avoid each other reflects the competition of ki- 
netic and potential energy. Such a topological (CO+BS) 
phase separation is believed to provide a minimization of 
the total energy as compared with its uniform supersolid 
counterpart. Thus, the parent checkerboard CO phase 
may gradually loose its stability under boson/hole dop- 
ing, while a novel topological self-organized texture is 
believed to become stable. The most probable possibil- 
ity is that every domain wall accumulates single boson, 
or boson hole. Then, the number of centers in a multi- 
center skyrmion nucleated with doping has to be equal 
to the number of bosons/holes. In such a case, we an- 
ticipate the near-linear dependence of the total BS vol- 
ume fraction on the doping. Generally speaking, one 
may assume scenario when the nucleation of a multi- 
center skyrmion occurs spontaneously with no doping. In 
such a case we should anticipate the existence of neutral 
multi-center skyrmion-like object with equal number of 
positively and negatively charged single skyrmions. How- 
ever, in practice, namely the boson/hole doping is likely 
to be a physically clear driving force for a nucleation of a 
single, or multi-center skyrmion-like self-organized collec- 
tive mode in the form of multi-center charged topological 
defect which may be (not strictly correctly) referred to 
as multi-skyrmion system akin in a quantum Hall fer- 
romagnetic state of a two-dimensional electron gas^S In 
such a case, we may characterize an individual skyrmion 
by its position (i.e., the center of skyrmionic texture), 
its size (i.e., the radius of domain wall), and the orien- 
tation of the in-plane components of pseudo-spin (U(l) 
degree of freedom) . An isolated skyrmion is described by 
the inhomogeneous distribution of the CO parameter, or 
staggered boson density l z , order parameter m z charac- 
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terizing the deviation from the half-filling, and m± that 
corresponds only to the modulus of the superfluid order 
parameter. 

It seems likely that for a light doping any doped parti- 
cle (boson/holes) results in a nucleation of a new single- 
skyrmion state, hence its density changes gradually with 
particle doping. Therefore, as long as the separation be- 
tween skyrmionic centers is sufficiently large so that the 
inter-skyrmion interaction is negligible, the energy of the 
system per particle remains almost constant. This means 
that the chemical potential of boson or hole remains un- 
changed with doping and hence apparently remains fixed. 

The multi-skyrmionic system in contrast with uniform 
ones can form the structures with inhomogeneous long- 
range ordering of the modulus of the superfluid order 
parameter accompanied by the unordered global phases 
of single skyrmions. Such a situation resembles in part 
that of granular superconductivity. 

In the long-wavelength limit the off-diagonal ordering 
can be described by an effective Hamiltonian in terms 
of U(l) (phase) degree of freedom associated with each 
skyrmion. Such a Hamiltonian contains a repulsive, long- 
range Coulomb part and a short-range contribution re- 
lated to the phase degree of freedom. The latter term can 
be written out in the standard for the XY model form 
of a so-called Josephson coupling 

Hj = - Jjj cos(ipi -ipj), (4) 

{id) 

where (pi, <pj are global phases for skyrmions centered at 
points respectively, Jij Josephson coupling param- 
eter. Namely the Josephson coupling gives rise to the 
long-range ordering of the phase of the superfluid order 
parameter in a multi-center skyrmion. Such a Hamil- 
tonian represents a starting point for the analysis of 
disordered superconductors, granular superconductivity, 
insulator-superconductor transition with array of 

superconducting islands with phases ipi,(pj. Calculating 
the phase-dependent part of skyrmion-skyrmion interac- 
tion Timm et alm& arrived at negative sign of Jij that 
favors antiparallel alignment of the U(l) pseudospins. In 
other words, two skyrmions are assumed to form a pe- 
culiar Josephson n micro-junction. There are a number 
of interesting implications that follow directly from this 
result : 2 i the spontaneous breaking of time-reversal sym- 
metry with non-zero supercurrents and magnetic fluxes 
in the ground state, long-time tails in the dynamics of the 
system, unconventional Aharonov-Bohm period hc/Ae, 
negative magnetoresistence. 

To account for Coulomb interaction and allow for 
quantum corrections we should introduce into effective 
Hamiltonian the charging energySi 

H ch = -^^riiiC-^ijUj , 

id 

where rij is a boson number operator for bosons bound 
in i-th skyrmion; it is quantum-mechanically conjugated 



to ip: m — —id/dipt, (C stands for the capacitance 
matrix, q for bosonic charge. 

Such a system appears to reveal a tremendously rich 
quantum-critical structure. In the absence of disorder, 
the T = phase diagram of the multi-skyrmion system 
implies either triangular, or square crystalline arrange- 
ments (Skyrmion crystal) with possible melting transi- 
tion to a Skyrmion liquid. It should be noted that anal- 
ogy with charged 2D Coulomb gas implies the Wigner 
crystallization of multi-center skyrmion with Wigner 
crystal to Wigner liquid melting transition, respectively. 
Naturally, that the additional degrees of freedom for 
skyrmion provide a richer physics of Skyrmion lattices. 
For a WC to be an insulator, disorder is required, which 
pins the WC and also causes the crystalline order to have 
a finite correlation length. Traditional approach to a 
Wigner crystallization implies the formation of a WC for 
densities lower than a critical density, when the Coulomb 
energy dominates over the kinetic energy. The effect of 
quantum fluctuations leads to a (quantum) melting of 
the solid at high densities, or at a critical lattice spac- 
ing. The critical properties of a two-dimensional lattice 
without any internal degree of freedom are successfully 
described applying the BKT theory to dislocations and 
disclinations of the lattice, and proceeds in two steps. 
The first implies the transition to a liquid-crystal phase 
with short-range translational order, the second does the 
transition to isotropic liquid. Disorder pins the Skyrmion 
lattice and also causes the crystalline order to have a fi- 
nite correlation length. For such a system provided the 
skyrmion positions fixed at all temperatures, the long- 
wave-length physics would be described by an antiferro- 
magnetic XY model with expectable BKT transition and 
gapless XY spin-wave mode. 

As regards the superfluid properties the skyrmionic 
liquid reveals unconventional behavior with two critical 
temperatures Tbs < t and T c < J, Tbs being the tem- 
perature of the ordering of the modulus and T c < Tbs 
that of the phase of order parameter 

The low temperature physics in Skyrmion crystals 
is governed by an interplay of two BKT transitions, 
for the U(l) phase and positional degrees of freedom, 
respectively. 2 * 3 Dislocations in most Skyrmion lattice 
types lead to a mismatch in the U(l) degree of freedom, 
which makes the dislocations bind fractional vortices and 
lead to a coupling of translational and phase excitations. 
Both BKT temperatures either coincide (square lattice) 
or the melting one is higher (triangular lattice) S 

Quantum fluctuations can substantially affect these re- 
sults. Quantum melting can destroy U(l) order at suf- 
ficiently low densities where the Josephson coupling be- 
comes exponentially small. Similar situation is expected 
to take place in the vicinity of structural transitions in 
Skyrmion crystal. With increasing the skyrmion density 
the quantum effects result in a significant lowering of the 
melting temperature as compared with classical square- 
root dependence. The resulting melting temperature can 
reveal an oscilating behavior as a function of particle den- 
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sity with zeros at the critical (magic) densities associated 
with structural phase transitions. 

In terms of our model, the positional order corresponds 
to an incommensurate charge density wave, while the 
U(l) order does to a superconductivity. In other words, 
we arrive at a subtle interplay between two orders. The 
superconducting state evolves from a charge order with 
Tc < T m , where T m is the temperature of a melting 
transition which could be termed as a temperature of 
the opening of the insulating gap (pseudo-gap!?). 

The normal modes of a dilute skyrmion system (multi- 
center skyrmion) include the pseudo-spin waves propa- 
gating in-between the skyrmions, the positional fluctu- 
ations, or phonon modes, of the skyrmions which are 
gapless in a pure system, but gapped when the lattice 
is pinned, and, finally, fluctuations in the skyrmionic in- 
plane orientation and size. The latter two types of fluc- 
tuation are intimitely connected, since the ^-component 
of spin and orientation are conjugate coordinates because 
of commutation relations of quantum angular momentum 
operators. So, rotating a skyrmion changes its size. The 
orientational fluctuations of the multi-skyrmion system 
are governed by the gapless XY model.— The relevant 
model description is most familiar as an effective theory 
of the Josephson junction array. An important feature of 
the model is that it displays a quantum-critical point. 

The low-energy collective excitations of skyrmion liq- 
uid includes an usual longitudinal acoustic phonon 
branch. The liquid crystal phases differ from the isotropic 
liquid in that they have massive topological excitations, 
i.e., the disclinations. One should note that the liquids 
do not support transverse modes, these could survive in 
a liquid state only as overdamped modes. So that it is 
reasonable to assume that solidification of the skyrmion 
lattice would be accompanied by a stabilization of trans- 
verse modes with its sharpening below melting transition. 
In other words an instability of transverse phonon modes 
signals the onset of melting. 

A generic property of the positionally ordered 
skyrmion configuration is the sliding mode which is usu- 
ally pinned by the disorder. The dcpinning of sliding 
mode(s) can be detected in a low-frequency and low- 
temperature optical response. 

In conclusion, the boson/hole doping of the hard-core 
boson system away from half-filling is assumed to be a 
driving force for a nucleation of a multi-center skyrmion- 
like self-organized collective mode that resembles a sys- 
tem of CO bubble domains with a Bose supcrfluid and 
extra bosons both confined in domain walls. Such a 



topological CO+BS phase separation, rather than an uni- 
form mixed CO+BS supersolid phase, is believed to de- 
scribe the evolution of hc-BH model away from half- 
filling. Starting from the classical model we predict the 
properties of the respective quantum system. In frames 
of our scenario we may anticipate for the hc-BH model 
the emergence of an inhomogeneous BS condensate for 
superhigh temperatures Ttps < t, and 3D superconduc- 
tivity for rather high temperatures T c < J < t. The 
system is believed to reveal many properties typical for 
granular superconductors, CDW materials, Wigner crys- 
tals, and multi-skyrmion system akin in a quantum Hall 
ferromagnetic state of a 2D electron gas. Topological in- 
homogeneity is believed to be a generic property of 2D 
hard-core boson systems away from half-filling. 

Despite all shortcomings, MFA and continuous approx- 
imation are expected to provide a physically clear semi- 
quantitative picture of rather complex transformations 
taking place in bare CO system with doping, and can be 
instructive as a starting point to analyze possible scenar- 
ios. First of all, the MFA analysis allow us to consider 
the antiphase domain wall in CO phase to be a very ef- 
ficient ring-shaped potential well for the localization of 
a single extra boson (hole) thus forming a novel type of 
a topological defect with a single-charged domain wall. 
Such a defect can be addressed as a charged skyrmion- 
like quasiparticlc which energy can be approximated by 
its classical value for CO bubble domain. It is of great 
importance to note that domain wall simultaneously rep- 
resents a ring-shaped reservoir for Bose superfluid. 

Unfortunately, we have no experience to deal with 
multi-center skyrmions as regards its structure, energet- 
ics, and stability. It should be noted that such a texture 
with strongly polarizable centers is believed to provide an 
effective screening of long-range boson-boson repulsion 
thus resulting in an additional self-stabilization. Nucle- 
ation of topological phase is likely to proceed in the way 
typical for the first order phase transitions. The present 
paper establishes only the framework for analyzing the 
subtleties of the phase separation in a lattice hc-BH 
model away from half-filling. Much work remains to be 
done both in a macroscopic and microscopic aproaches. 
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